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Abstract 



We find a representation of smooth solutions to the Cauchy problem for a scalar 
multidimensional conservation law as small diffusion limit of a stochastic perturba- 
tion along characteristics. It helps, in particular, to study the process of singularities 
formation. Further, we introduce an associated system of balance laws that can be 
interpreted as describing the motion of a continuum with some specific pressure 
term. This term arises only after the instant when the solution to the initial Cauchy 
problem looses its smoothness. Before this instant the system coincides partly with 
the one known as pressure free gas dynamics. 
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Introduction 

We consider the initial value problem 

n 

u t + ^(h(t,x,u)u Xi = Q, u{0,x) = u {x), u (x) G Cl(R n ;R), (1) 

i=l 

where t G M+, x G M n , di(t,x,u), i = l,...,n, is a real-valued C l function 
defined on some open subset of (R + x M™ x M). For a technical reason the 
functions di(t, x, u) are assumed to grow at infinity not quicker than a linear 
function in u. 

A particular important special case is given by the scalar conservation law in 
the form 

u t + divF(t, u) = 0, 

where F(t,-) = (Fi(t, ■),..., F n (t, ■)) is a C 2 vector-function defined on some 

du ' 



open subset of E, for any t G R+, and aj(t, tt) = — diw^ ^ = 1, ..., n. 



The main aim of this paper is to obtain an asymptotic formula for the solution 
of the Cauchy problem ([1]) for the case of a scalar conservation law. The for- 
mula is obtained by the limit for vanishing perturbation of the corresponding 
stochastically modified equation (small diffusion limit). 

Nevertheless, let us first consider the general case. 

Let us write the associated characteristic ODE: 

dxi du . 

— — = aAt.x.u), -— = 0, t = l,...,n. 
dt K ' h dt 

Its stochastic analog is 

dXi{t) = ai {t, X(t), U{t))dt + (TtdiWDt, dU(t) = a 2 d(W 2 ) t , (2) 

Xi(Q) = x { , U{0) =u, t > 0, 

i = l,...,n, X(t) and U(t) are considered as random variables with given 
initial distributions, (X(t), U(t)) runs in the phase space M n x]R 1 , o\ and a 2 are 
nonnegative constants such that |cr| ^ (a = (ai,a 2 )) and ((V^ 1 )^, (W 2 ) t ) = 
{Wi, . . . , W^, W 2 )t is an n+1 - dimensional Brownian motion, i.e. the Wl, W 2 , 
i = 1, ...,n, are independent one-dimensional standard Brownian motions. 

Let P(t,dx,du), t G M + , x G M n , be the probability of the joint distribution 
of the random variables (X, U), subject to the initial data 

Po(dx, du) — S u (uo(x)) po(x)dx, (3) 



where po is a bounded nonnegative function from C(M n ) and dx is Lebesgue 
measure on M. n , 5 U is Dirac measure concentrated on u. P(t, dx, du) has the 
form P(t, x, du)dx, where P(t, x, du) is a positive measure with respect to u 
and a function with respect to x (density function of P(t, dx, du) with respect 
to Lebesgue measure). 

We look at P — P(t, dx, du) as a generalized function (distribution) with 
respect to the variable u. It satisfies the Fokker-Planck equation 



dP 

~dt 



n Q n \ d 2 n \ d 2 

E T^T a ^' x > u ) + E o a i 2 7^2 + E o a 2 



fc=l aXfc fc=l A OX k k=l Z ° n fe 



P, (4) 



subject to the initial data (j3j) 



There is a standard procedure for finding the fundamental solution for (jlj) (see, 
e.g. [9]). This procedure consists in a reduction of the equation to a Fredholm 
integral equation, the solution of which can be found in the form of series. 
We are going to show that for a(t, x, u) = a(t, u) one can also find an explicit 
solution to the Cauchy problem (JH), ©• 



Let us introduce, still in the general case, the functions, for t G R+, x G 
depending on a = (cr^, cr 2 ): 



r, 



p ff {t,x) = / P(t,x,du), (5) 



/ uP(t, x, du) 
U ^ x) = R fP(t,x,du)> (6) 



J a(t,x, u)P(t, x, du) 
a - (t ' x) = R JP(t,x,du) ' (?) 

the integrals in the numerator being assumed to exist in the Lebesgue sense. 

It will readily be observed that u„(0,x) = uq{x) and a a (0,x) = a(0,x,Uo(x)). 

We denote 



p(t,x) = lira p a (t,x), u(t, x) = lim u a {t, x), a(t, x) — lim a a (t,x), 

a— >0 a— >0 cr— >0 

provided these limits exist. 



1 Case of a conservation law 



Now we dwell on the simpler case of a conservation law, where a = a(t, u). Here 
the equation (J3J) can be solved explicitly. Moreover, for the sake of simplicity 
we set a 2 = and denote o\= o. 

Proposition 1 If a = a(t,u), then problem (J3J), ([3]) has the following solution: 

n t 

y] \J a i( T , u o(y)) dT +yi- x i\ 

If _i = l 

P(t,x,du) = ■== / S u (uo(y))po(y)e ^ dy, (8) 

V 2%ta) n J 

t > 0, x G W" 1 , or, in other words, 

n t 

/] J a i( T ,uo(y)) dT + Vi- x i\ 

1 r »_i o 

<f>(u) P(t, x, du) = -r-7=-r- / (p(u (y)) p (y) e ^i dy, 

V Into •)" J 

(9) 
for all (^(u) eC (M). 

Proof. We act as in [3], [I]. Namely, we apply the Fourier transform to 
P(t,x,du) in (J3J), (j3j) with respect to the variable x and obtain the Cauchy 
problem for the Fourier transform P = P(t, X, du) of P(£, x, du): 

r)P 1 

^ = -(-a 2 |A| 2 + z(A,a(^)))P, (10) 

P(0,X,du) = [ e- i ^8 u (uo(y)) Po (y)dy, XeR n . (11) 



Equation ffTUl) can easily be integrated and we obtain the solution given by 
the following formula: 

t 

-±a 2 \\\ 2 t+i f(\,a(T,u))dr 

P(t,X,du) = P(0,X,du)e o . (12) 

The inverse Fourier transform (in the distributional sense) allows to find the 
density function P(t, x, du), t > 0: 

Pit, x, du) = —V / e i(x ' x) P{t, A, du) d\ = 
(27n n J 



(2n] 



2n 



J(X,x) 



( 



t 



e 



—i f(X,a(r,u))dT 

^e o 5 u (u (y)) Po (y)dy\ e^^dX 



\oH 



(27T) 



-T- / 8 u (My)) Po(y) / e V 



/ r A' , r 

i\x — I a(T,u)dT — y\ \ \ I a(r y u)dr 


~~7*t 

/ 



A-- 



+ v — x[ 



dXdy 



1 



2Trta) 



\\ a(T,UQ(y))dT + y~x\ 

5 u (u (y)) p (y) e~~ ** dy, t>0,xeR n . 



The third equality is satisfied by Fubini's theorem, which can be applied by 
the absolute integrability and the bound on the function involved. Thus, the 
proposition is proved. 

Remark 1 In the general case o<i ^ an analogous formula can be obtained 
in a similar way. 

Corollary 1 The functions p a , u a and a a defined in (jSJ) - (ED can be repre- 
sented by the following formulae: 



/] J a i( T ,' u o(y))dT + y i -x i \ 
i=l 

pa(t,x) = / p (y)e 5^ dy, 



(13) 



u a (t,x) 



n t 

}, I J a i (T,u (y))dr + y i -x i \ 
i=l 

Ju (y)p Q (y)e &* dy 

R« 

n t > 

^ I J a i(- r ,u (y))dT + y i -x i \ 2 

/ Po(y) e ^ dy 



(14) 



U//T I Ul ti-l J 



/] \J a i(i~, u o(y)) dT +yi- x i\ 2 

i-l 

/ a(t,uo(y))po(y)e ^ ds 



I Po(y) e 



n t 

y^ i j a i( T .«o(«)) dT 

8 = 1 



'ds 



(15) 



Proof. The result is obtained by substitution of P(t,x,du) as given by 
in ©, © and Q. 



1.1 Asymptotic formula for smooth solutions 
Let us define the following subset A of R: 

i n 

teA if yeiJ H^)«( r ' M °^))M^A > - 1 : (16) 



i=1 



where u G C^(R n ). It is not difficult to show the A is an open set. Denote 
t*(u ) = sup A. 

The following theorem holds: 

Theorem 1 Let u(t, x) be a solution to the Cauchy problem 

n 

u t + ^2ai(t,u)u Xi = 0, u(0,x) =u (x), (17) 

i=i 

where a{, i = 1, ...,n, are C 1 - functions defined on some open subset of (R + x 
R) and u G C^(R n ). Assume £*(uo) = sup A > 0, A 6em^ defined by n~b}) . 
Then for t G [0, £*(ii )), 

ti(t, x) = -u(t, x) = lim u a (t, x), 

a— >0 

where u a (t,x) is given by (0|) and the limit exists pointwise. 

Proof. The proof is similar to the one given in [1] for a related problem. 
According to the classical theory (see, e.g. [6]), Theorem 5.1.1), the solution 
u of (IT7|) exists on some maximal interval [0, T), T < oo and is a C l - smooth 
function. Since u is constant along characteristics, its value at any point (t, x), 
with x G R n , t G R+, satisfies the implicit relation 

t 
u(x,t) = u (x — / a(T,u)dr). (18) 

o 

In particular, the range of u coincides with the range of Uq. 
Differentiating (1181) yields 

d yi u (y) 



d Xi u(t,x) = — , y = x- a(r,u)dr. (19) 

1 + / T,((k) u (r,uo(y))(uo(y)) Vi dT o 

«=1 

This imply T = t*(uo). If < £*(w ) < +oo, then the solution to the Cauchy 
problem blows up at the instant t*(uo). Otherwise, the solution keeps its 
smoothness for all t > 0. 



The formula (Q implies, using the weak convergence of measures and the fact 
that p and u are continuous and bounded and independent of a 



J a("<~,uo(y))dT + y- 



II 



lim u ff (t, x) 



I u (y) P o(y) lim jj^e- =* dy 



cr^O " v ' ' * 



I / a(T,u (y))dT+y-x\ 2 







l*MHv£zF e ^ dy 



I uo(y)po(y)8 p (t,x,y)dy 



I Po(y)$ P (t,x,y)dy 
with 

^ y) :-l^^r + y -, ( 2 0) 



where 5 P is the Dirac measure at p G R n . We can use locally the implicit func- 
tion theorem and find y = yt, x (p) from p(t,x,y). The condition for existence 
of this function is the invertibility of the matrix 

Cij{t,y) = — l -^ — , i,j = l,...,n. 

9yj 

This matrix fails to be invertible for t = t* (u Q ) . For t < t* (ito) 

u(t, x) = lim u ff (t,x) = 

m [ Uo(yt,x(p))po(ytM) det(C(t,y 4 , x (p)))" 1 5 p (dyt, x ) 

— = u (yt,x(°))- 



I Po(yt,x(p)) det(C(t,|/ ti:r (p)))- 1 5 p (dy tiX ) 



Let us introduce the new notation y Q (t,x) = y tiX (0)- Then f )20l) implies the 
following vectorial equation: 

a(T,Uo(y (T,x)))dT + yo(t,x)-x = 0, t > 0, iGf". (21) 



Let us show that u{t,x) = u (y (t,x)) satisfies equation ([T]), that is 

n n 

Yj d j( u o)(yo,j)t + J2 ai(*>«o)3fe(«o)(2/o,k)ay =0. (22) 

j=l j,k=l 

and Wo (2/o(0, x)) = Wo(^)- Here we denote by y ,j the z - th components of the 
vector y . 



For t < t*(«o) we can differentiate ([21]) with respect to t and Xj to get the 
matrix equations: 

n 

J2 C v (Vo,j)t + %,i = 0, i = 1, ..., n, 

i=i 

and 

n 

XI C*k (yo,k)xj + Sij = 0, i,j = 1, ...,n, 
fe=i 

where <5jj is the Kronecker symbol. The equations imply 

n 

(yo,j)t = - X)(C _1 )ijUo,i, (j/o,fc)^ = -(C -1 )jfc- (23) 

8=1 

It remains now only to substitute (1231) into (l22j) to see that w(£, x) satisfies the 
first equation in ([IT]) . 

Further, (J2"T|) implies m (?/o(0,x)) = %(£), thus Theorem [1] is proved. 

Remark 2 For a^ = aj(w) (i.e. cii is independent of variable t) we have the 
Conway' criterium $5$: 



t*(u ) = sup 



/ \ 

1 

n 

y E(ai)u(uo(y))(uo(y)) yi J 

Note that if the denominator vanishes, then t*(ito) = oo and the solution does 
not blow up. If t*(uo) < 0, then the solution is globally smooth for t > 0, as 
well. 

Proposition 2 Under the assumptions of Theorem^ the vector 

a(t, x) = lim a a (t, x), 

a— >0 

where a a (t,x) is given by |7|) ; solves the multidimensional Burgers equation 

(a(t, u))t + (a(t, u), V)a(t, u) — 

with initial data d(0,u(0,x)) = a(0,uo(x)). 

Proof. This fact follows directly from Proposition 2.1 of [I]. 

Remark 3 The introduction of a small perturbation of deterministic equation 
to study then the original equation in the limit of vanishing noise has appeared 
in several contexts, particularly for equations of the reaction- diffusion type, 
see, e.g. JT$, J?]/ and references therein. 



1.2 Associated system of balance laws 



Now we consider the following question: what system of equations do the 
triple (p ai u a) a a ) and its limit (p,u,a) satisfy before and after the blow up 
time i*(«o)? 

The following proposition holds: 

Proposition 3 The functions p a , u a anda a , given by (J5J), flOJ) and (j7|), satisfy 
for t > the following PDE system: 

£ + **■<**> = k£7*r (24) 

^> + div.fe «. <v) = ^» £ %p> - 5, (25) 

(u - M CT (t, a?))((a(t, w) - a„(t, x)), V x P{t, x, du)); 



<9(p CT a ffii ) 1 ™ <9 2 (p CT a CT ,;) 

— + dw x {p a a^a a ) = -a ^ — -^ J^, (26) 

i = 1, ..,n, where 

I° { = (ai(t,u) - a^i(t,x)((a(t,u) - a rT (t,x)),V x P(t,x,du))+ 



ai(t,u)) t P(t,x,du). 



Proof. The equation (1211) follows from the Fokker-Planck equation (j3j) di- 
rectly. 

Let us prove (I2"6"j) (the derivation of (1231) is analogous). We note that the 
definitions of a a (t,x) and p a (t,x) imply 



dip an) d r f 

— — ■ — = — / a(t,u)P(t,x,du) = / a(t,u)P t (t,x,du) = 

- J a(t, u)(a(t, u), V x P(t, x, du)) + -a 2 £ ^^, (27) 



fc=i ^fc 



where P t = ^P. 



Further, we have 

a a ' k a ' 1 = a ffli (t, x) - — / a k (t, u) P(t, x, du)) + 

ox k dx k \ R n j 

r\ I j ai(t,u) P(t,x,du) 

/ a k (t, u) P(t, x, du) — R " j p(t|3 . |du) — | = I a a ,i(t, x) a k (t, u) P Xk (t, x, du) 4 



I ai(t,u) P x ,(t,x,du) I P(t,x,du) — I a,i(t,u) P(t,x,du) I P x ,(t,x,du) 

J a k {t,u)P(t,x,du)^ ^ 



J P(t,x,du) ) 



/ (a k (t,u)a a)i (t,x) + ai(t,u)a a , k (t,x) - a a . k (t,x) a ff) i(t,x)) P Xk (t,x,du), 

i,k = l,...,n, with P Xh = -^P. 

Equation (126|) follows immediately from (1271) and (|28|) . Thus, Proposition [3] is 
proved. 



Corollary 2 Before the instant t*(uo), the blow up time of the solution to the 
Cauchy problem ( TT71) . the triple (p, u, a), which constitutes the limit as \a\ — > 
of the triple (p CT , u a , a a ), solves the following system: 

d t p + dw x {pa) = 0, (28) 

d t {pu) + V x {pua) = 0, (29) 

dt(pa)+V x (pa®a) = Q. (30) 

Proof. Equation ( |28l) follows from the properties of parabolic differential 
equations with a small parameter in front of the derivatives of second order ( [S j , 
Theorem 3.1), since until the instance t*(uo) the coefficients of equation (|24j) 
are differential) le. Equation (|29l) follows from (1281) and Theorem [TJ Proposition 
M implies ((30]). 

Remark 4 System (|28|) and (1301) constitutes the so called pressureless gas 
dynamics system , the simplest model introduced to describe the formation of 
large structures in the Universe, see, e.g. [10]. 

Remark 5 As it has been shown in OV on an example, for discontinuous 
solutions to (TT71) the limits as a — > of the terms J" and 1% do not vanish as 
a — > and yield some specific pressure. 



10 



Remark 6 The method of special stochastic perturbations, applied here, was 
used in fBjj, |5j/ for studying other deterministic problems. 



References 



[1] S.Albeverio, L. Di Persio, E. Mastrogiacomo, Small noise asymptotic 
expansions for stochastic PDE's, I. The case of a dissipative polynomially 
bounded non linearity, Tohoku Mathematical Journal, 63(2011), 877-898. 

[2] S. Albeverio, O. Rozanova, The non-viscous Burgers equation associated with 
random position in coordinate space: a threshold for blow up behaviour, 
Mathematical Models and Methods in Applied Sciences, 19(2009), No. 5, 1-19. 

[3] S. Albeverio, O. Rozanova, Suppression of unbounded gradients in SDE 
associated with the Burgers equation, Proceedings of Amer.Math.Soc. 
138(2010), 241-251. 

[4] S. Albeverio, A. Korshunova, O. Rozanova, Probabilistic model associated with 
the pressureless gas dynamics, (2009), arXiv:0908.2084, submitted. 

[5] E.Conway, The formation and decay of shocks of a conservation law in several 
dimensions, Arch. Rat. Mech. Anal, 64(1977), 47-57. 

[6] C.Dafermos, Hyperbolic conservation laws in continuum physics, Springer, 
2010. 

[7] F.Flandoli, Random perturbation of PDEs and fluid dynamic models. Lectures 
from the 40th Probability Summer School held in Saint-Flour, 2010. Lecture 
Notes in Mathematics, 2015. Springer, Heidelberg, 2011. 

[8] M.I.Freidlin, A.D.Wentzel, Random Perturbations of Dynamical Systems, 
2nd ed., Springer, New York, 1998. 

[9] A. Friedman, Partial Differential Equations of Parabolic Type, R. E. Krieger 
Publishing Company, Malabar, Florida, 1983. 

[10] S.Shandarin, Ya.B.Zeldovich, The large structure of the universe: turbulence, 
intermittence, structures in a self-gravitating medium, Rev. Modern. Phys, 
61(1989), pp. 185-220. 



11 



